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Abstract. The ABn-i chain is a system that consists of repeating a unit cell with N sites, where between 
the A and B sites there is an energy difference of A. We considered bosons in these special lattices and 
took into account the kinetic energy, the local two-body interaction, and the inhomogenous local energy 
in the Hamiltonian. We found the charge density wave (CDW) and superffuid and Mott insulator phases, 
and constructed the phase diagram for N = 2 and 3 at the thermodynamic limit. The system exhibited 
insulator phases for densities p = a/N, with a being an integer. We obtained that superffuid regions 
separate the insulator phases for densities larger than one. For any N value, we found that for integer 
densities p, the system shows p -I- 1 insulator phases, a Mott insulator phase, and p CDW phases. For 
non-integer densities larger than one, several CDW phases appear. 

PACS. 05.30.Jp Boson systems - 05.30.Rt Quantum phase transitions 


1 Introduction 

The study of the ground state of bosonic systems is an 
interesting area within the study of cold atoms due to the 
capability of emulating them in optical lattices, with the 
advantage of having absolute control over the parameters 
(kinetic energy, interactions, and density) and access to 
various dimensions m- 

Experimental progress on this topic has allowed the 
prediction of the superfluid-Mott insulator transition in 
ultracold bosonic atoms [3] and the observation of quan¬ 
tum phase transitions [4] in an optical lattice, which has 
opened the way for potential studies that have revealed 
the fundamental physics of these systems. 

Recently, the implementation of optical lattices has 
motivated the exploration of bosonic systems in super¬ 
lattices whose arrangement is characterized by a periodic 
potential [5116] . for instance an energy difference of A in 
the unit cell between sites called A and TV — 1 B in the 
unit cell (notation ABf^-i). Bidimensional systems with 
AB configurations have been created confining ®^Rb [7] 
and m atoms in optical lattices. This experimental 
progress allows us to believe that experimental study of 
the one-dimensional ABjv-i chains will be done in the 
coming years. 

In 2004, Buonsante and Vezzani elaborated an an¬ 
alytical description of the physics of cold bosonic atoms 
trapped in super lattices. Within the study of the finite- 
temperature phase diagram, they included the results for 
zero temperature, and insulator domains for fractional 
densities were demonstrated. The phase diagrams due to 
the interplay between the on-site repulsive interaction, su¬ 



Fig. 1. (color online). Schematic representation of the super¬ 
lattices considered in this paper. The chain AB (a) is formed 
by the superposition of two waves, where the first has double 
the frequency of the other, while the chain ^152 in (b) is formed 
by the superposition of two waves where the first has triple the 
frequency of the other. We consider that the difference of the 
hopping parameter between neighboring states is very small. 


perlattice potential strength, and filling were found by 
Rousseau et al. m- 

Later, Dhar and colleagues [II1II1II3] studied the quan¬ 
tum phases that emerge when the AB chain is considered 
(see Fig. [IKa)). Such studies proved that a new insulator 
phase results at density p = 1/2, and the Mott insulator 
phase corresponding to p = 1 undergoes a phase transi¬ 
tion to another insulator when the parameter A is near 
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Fig. 2. Graphs of chemical potential versus 1/L for X/t = 20 
in a) p = 1 and b) p = 3/2. The points are DMRG results, and 
the lines represent the regression to the thermodynamic limit. 


the strength of the interaction between the particles (t/), 
leading to a different boson arrangement. 

Although the results of these new phases indicate the 
consequences produced by the superlattice on a boson 
chain, we believe that there are new aspects to be dis¬ 
covered. For instance, we can ask: is it possible that the 
insulator phases at p > 1 evolve to other insulator phases 
passing briefly through a superfluid when A increases? 
What happens in the phase diagram if we increase the 
number of sites in the unit cell? To answer these ques¬ 
tions, we used the density matrix renormalization group 
method m and the von Neumann entropy m to study 
the quantum phases of the model. We found that for frac¬ 
tional densities lower than one there is a unique CDW 
phase for any value of A. For densities larger than one, 
we obtained a finite number of insulating phases when A 
increases, two CDW phases for fractional densities, and 
p -|- 1 for integer ones. The von Neumann entropy shows 
that the insulator phases are separated by superfluid re¬ 
gions. 

This paper is organized in the following way: in sec¬ 
tion 2, we explain the Hamiltonian for AHjv-i chains. The 
thermodynamic limit phase diagrams and the von Neu¬ 
mann entropy results for AB and AB 2 chains are shown 
in sections 3 and 4. Finally, the conclusions are presented 
in section 5. 


2 Model 

The ground state of ultracold bosonic atoms in homoge¬ 
neous systems is described by the Bose-Hubbard model. 
When a system of bosons in an inhomogeneous lattice is 
considered, the Hamiltonian is given by: 

H = {a\aj + H.c.) + ^ X! “ 1) + X! 
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Fig. 3. (color online). Phases and transitions in the AB chain. 


( 1 ) 

where t is the hopping parameter, < i,j > denotes a pair 
of nearest-neighbor sites i, and j, aj (a^) creates (annihi¬ 
lates) a boson at site i. U represents the local interaction 
in the second term of the Hamiltonian, where ni = aj ai is 
the number operator, and Ai denotes the shift in the en¬ 
ergy levels of the sites in each unit cell. We set our energy 
scale taking t = \ and the interaction parameter U/t = IQ. 

Particle chains in a super lattice are denoted as AB^-i, 
where the site A has a difference of energy of A from the 
iV — 1 sites called B per unit cell. In the present paper, 
we consider two types of superlattices. The first superlat¬ 
tice has a potential with periodicity 2 {AB chains), mean¬ 
ing two sites per unit cell with a potential difference A, 
schematically represented in Fig.lIJa). The second one has 
a periodicity of 3 and is shown in Fig. [T](b) {AB 2 chains). 

As is well known, in a homogeneous environment the 
one-dimensional bosonic system exhibits a phase transi¬ 
tion between a Mott insulator phase, characterized by an 
integer filling, and the superfluid phase, which is compress¬ 
ible [TFlITTlfTB] . In a superlattice type AB, the behavior of 
the system in the ground state exhibits additional phases 
for half-filling and density equal to one, where insulating 
phases due to the superlattice potential are exhibited m- 
The current investigation involved calculating the en¬ 
ergies E{N, L) for lattices with different lengths L and N, 
-I- I and — 1 particles such that we obtain the chem¬ 
ical potential for the increase (/ip) and decrease (/i/j) the 
number of particles to one, where the general expressions 
are: 

fip = E{N + 1,L)-E{N,L) (2) 

and 

Hh=E{N,L)-E{N-l,L). (3) 

To obtain the energies, we used the density matrix renor¬ 
malization group method for lattices from 24 to 84 sites, 
obtaining an error of around Then we extrapolated 

our results at the thermodynamic limit and repeated the 
process for a wide range of values of the parameter A. For 
example, all insulator phases are characterized by a gap. 
The analysis of this quantity at the limit L —>■ 00 gives us 
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Fig. 4. (color online). Density profile p versus chemical poten¬ 
tial p for bosons in superlattice type AB and various values of 
A, with U/t = 10. 


Fig. 5. (color online). On-site number density plotted against 
lattice site index with L = 50 sites at density p = 3/2 and 
p = 2. The lines are visual guides. 


information about the phase in which the system exists. 
The gap for any boson system is given by: 

A = pp-ph- (4) 

If Z\ 0, the system is in an insulator phase; if it is not, 
the system is in a superfluid phase. This is shown in Fig. 
[2] for the AB chain for two cases: in Fig.[2Ka), we consider 
a density p = 1 and X/t = 20 and we observe that the sys¬ 
tem has a large gap at the thermodynamic limit, which 
indicates that the system is in a CDW phase, and specifi¬ 
cally the charge distribution generated by the superlattice 
structure is {2,0, 2,0, 2,0...}, as reported before (see Fig. 
0). If we maintain the A parameter the same and change 
the global density to p = 3/2, we expect that the system 
will be in a CDW state due to the superlattice, according 
to the results of Buonsante and Vezzani[^, with a den¬ 
sity profile {2,1, 2,1, 2, l...}(see Fig. 0. However, in Fig. 
0b), we observe a gap equal to zero at the thermodynamic 
limit, which indicates that the system is in a superfluid 
state. In conclusion, for X/t = 2U, the ground state is 
CDW for p = 1 and is superfluid for p = 3/2. 


3 AB phase diagram 

Dhar et al. m have shown that considering bosons in a 
superlattice with periodicity equal to two, a CDW phase 
with one particle per unit cell appears for p = 1/2 (see 
Fig. 0, which remains for any value of A, while for p = 1, 
they obtained a Mott insulator phase and a CDW phase, 
because when the hopping and the local interaction pa¬ 
rameters are fixed, for values oi X < U the ground state 
will be a Mott insulator with one particle per site. How¬ 
ever, for values of A ~ H, the local repulsion can be com¬ 
pensated by the superlattice term and the quantum fluc¬ 
tuations increase, delocalizing the particles, which leads 
to a superfluid state. As A is increased further, the bosons 
prefer to be in the wells and the ground state is a CDW 


(see Fig. 0. They showed that the superfluid region be¬ 
tween the Mott insulator and the CDW phases is in the 
range {9.6,10.0}. In this section, we extend these results 
by considering densities larger than p = 1 in the AB chain. 

In Fig. 0 we show the density p versus the chemical 
potential p/t for some specific values of X/t, considering 
U/t = 10, because the previous phase diagram was con¬ 
structed for this value. The case of a homogeneous lattice, 
with X/t — 0, is displayed in Fig. 0a). In this figure, it 
is possible to clearly observe plateaus at integer densities 
of p = 1 and p = 2. The width of the plateaus indicates 
the size of the energy gap at the thermodynamic limit for 
each density. Between these values of p/t, a Mott insula¬ 
tor phase is found. For different values of p/t, the system 
is compressible and is in a superfluid phase. 

When an inhomogeneous lattice is considered, for ex¬ 
ample with A = 6, shown in Fig. 0b), the number of 
plateaus increases. It is crucial to emphasize that the plateaus 
for p = 1 and p = 2 are smaller than for A = 0, where 
the plateau for p = I is larger than the plateau for p = 2. 
Also, the other plateaus for half-integer p are shown as pre¬ 
dicted. These plateaus are originated by the superlattice 
structure type that characterizes the boson chain, while 
the plateaus for integer densities are caused by the inter¬ 
action between the particles of each site. It is important to 
highlight the fact that between each plateau at a specific 
value of X/t, we always obtain a superfluid region. 

In Fig. 0c), for X/t = 10, plateaus for densities p = 
1/2 and p = 3/2 are observed, while for p = I and p = 2 
the width of the plateaus is very small, and around A ~ 
we found that the Mott insulator phase disappears for 
all integer densities and a superfluid region anticipates a 
new insulator phase. On the basis of Fig. 0 for p = 2, 
we can explain why the Mott insulator phase disappears. 

A Mott insulator phase means p particles per site, which 
implies an energy Up{p — 1) + pA for two sites. For A « 

U, this energy can be compensated by a particle jump 
into wells, the energy for two sites is 17(p^ — p + 1) -b 
(p — I)A, and the fluctuations lead to a superfluid state 
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Fig. 6. (color online). Phase diagram of one-dimensional 
bosonic system for the AB chain. The points (open circles) 
represent the boundaries and are DMRG results. 


Fig. 7. (color online). Density versus the chemical potential 
for some values of \/t, for a system of bosons in a superlattice 
type AB 2 ■ 


for any density p. Note that the hopping term is small 
{U/t = 10). This contributes to the fluctuations, but is 
not the most important factor. This fact should be taken 
into consideration throughout this paper. 

For X/t = 16, the plateaus reappear, and we get insu¬ 
lator phases at all commensurate densities p= 1/2,1,3/2 
and p = 2. Also, we can observe that the plateau for 
p = 1/2 always increases as X/t increases, which is shown 
in Fig. SKd). In Fig. SKe), for X/t = 20 the plateau at 
p = 3/2 is not present, and the insulator phase due to the 
superlattice disappears around A ~ 2U. From Fig. |31 we 
see that the CDW phase for this density consists of two 
bosons in the wells and one in the barriers. Note that if 
the barrier particle jumps to the well, the energy for two 
sites 3C/ compensates for the energy of the CDW arrange¬ 
ment U + X, when A ~ 2U, which leads the system to a 
superfluid state. 

Finally, for X/t = 24, it is possible to observe that 
again we have plateaus for all the densities, and a new 
insulator phase for p = 3/2 is obtained, as is shown in 
Fig.SKf). 

The above discussion allows us to conclude that for 
all densities p > 1, a new CDW insulator phase occurs for 
larger values of X/t. For integer densities, this CDW phase 
is achieved by first passing through a Mott insulator phase 
and a superfluid one, while for p = 3/2, we obtained two 
CDW phases separated by a superfluid phase. 

Each of the above insulator phases that have been 
mentioned is characterized by a particular charge distri¬ 
bution in the system. In Fig. [SI the lattice density profile 
is shown for densities p = 3/2 and p = 2, with some pa¬ 
rameters taken before and after the transition region. In 
Fig. ISKa), for p = 3/2 and X/t = 10, the distribution of 
two bosons at one site and one boson at the other site is 
observed, while for X/t = 24, in Fig. |S{b), the charge dis¬ 
tribution is close to {3, 0,3, 0,3, 0...} . Clearly, the system 
experiences a change of state caused by the spatial struc¬ 
ture of the lattice. We emphasize that this possibility was 
not observed by Dhar et al. 


In Fig. I^c) for X/t = 2, a Mott insulator phase with 
two particles at each lattice site is shown, but as X/t in¬ 
creases to X/t = 24 in Fig. liKd), the charge distribution 
in the system indicates that the system is organized as 
{3,1,3,1,3,1...}. 

As was reported, the insulator phase at density p = 1/2 
for A/t yf 0 is always present, and for p = 1, a phase tran¬ 
sition is located near X/t ss U (see Fig. [5)). Also, this figure 
shows that for all commensurate densities with p > 1, the 
system exhibits a superfluid region between two insula¬ 
tor regions, and the parameter A can generate a quan¬ 
tum transition between these states due to the superlat¬ 
tice structure. Within the values of X/t considered, we see 
that there are two insulator phases for densities p > 1, and 
the gap of the CDW phase for the densities p = 1/2 and 
p = 1 saturates. Because for these densities, the barriers 
are unoccupied (see Fig.jH]), and this implies that no state 
change may occur. 


4 AB 2 phase diagram 

Now the case of bosons in a superlattice type AB 2 will 
be studied. In this superlattice, each unit cell is formed 
by three sites. Two i? sites differ by a potential of A with 
respect to A (see Fig. [T}b)). Note that real systems with 
this structure in one or three dimensions have been created 
and studied [Biiniiii]- 

Again, we consider a local repulsion interaction be¬ 
tween the bosons oi U/t = 10, and the chemical potential 
is determined at the thermodynamic limit. 

The evolution of the plot of density versus chemical 
potential as a function of the superlattice parameter is 
shown in Fig. [71 For X/t = 6 (Fig. |Z}a)), the number of 
plateaus increases, in contrast with the superlattice type 
AB. The first two plateaus are related to the possibility of 
having one or two bosons in the wells of each unit cell (see 
Fig.jS]), which means CDW phases with densities p = 1/3 
and p = 2/3. As in the AB case, we expected that these 
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Fig. 8. (color online). Phases and transitions in the AB 2 chain. 


plateaus would increase with A until they saturate. Also, 
we obtained other plateaus at commensurate densities of 
multiples of 1/3, as was predicted. 

Fig. [T^b) shows the results for X/t = 9.1, where the 
plateaus atp=l,p = 4/3 and p = 2 disappear, while 
the size of the gaps for the other commensurate densities 
remains finite. The reappearance of these plateaus hap¬ 
pens for larger values of A, as can be seen in Fig. EKc). 
This means that there is a superfluid region separating 
two insulator regions for A around U for each density. For 
integer densities, the explanation is similar to the previous 
case {AB); fixing the density p (each site in the well has p 
bosons) and a barrier strength A, the energy of two sites 
(the barrier and one of the well) is Up{p—l) + pX (see Fig. 
0). If one boson jumps into the well, previous energy can 
be compensated when X x U. Now the energy of the two 
sites is U{p'^ — p + 1) + {p — 1)X. We can conclude that for 
any AB]\[ chain, there is always a superfluid region around 
U separating a Mott insulator phase and a CDW one for 
any density p. 

Finally, in Fig. (T^d), with X/t = 19, the calculations 
show that the gaps for p = 5/3 and p = 2 are very small, 
which indicates that for these densities a superfluid region 
will appear around 2U. Note that in the AB chain, a su¬ 
perfluid region around 2U appears for the commensurate 
density p = 3/2, but this does not occur for the integer 
density p = 2 (see Fig. [6]). For X = 2U and p = 2, the 
charge distribution of the ground state is three particles 
in the well and one in the barrier (see Fig. 0, so the quan¬ 
tum transition happens when the particle in the barrier 
jumps to the well. The above happens for A Ri 3U. Before 
the transition, the energy for two sites is 3C/ -I- A, which 
can be compensated with the energy 6U once the particle 
has jumped into the well. The location of this transition 
marks one of our new hndings due to the increase of the 
number of sites of the wells. 

In order to distinguish between insulators, we calcu¬ 
late the arrangement of the particles in the superlattice 
according to the increase of the parameter X/t; i. e. the 
density profiles for values oi X/t lower and higher than the 



Fig. 9. (color online). On-site number density plotted against 
lattice site index at densities p = 4/3, 5/3 and p ~ 2 and values 
smaller or larger than the critical point. The lines are visual 
guides. 


transition region are shown in Fig. [HI where the densities 
p = 4/3, p = 5/3 and p = 2 are considered. For den¬ 
sity p = 4/3 in Fig. [Hja), with X/t = 2, the superlattice 
induces fluctuations of charge such that there is one par¬ 
ticle in the barrier and three in the well (see Fig. 0, while 
Fig. 0b), with X/t = 15, shows the particles organized 
such that the occupancy is {2, 2,0, 2, 2, 0,...}. The quan¬ 
tum phase transition takes place around U, which can be 
explained by taking into account that the smallest num¬ 
ber of particles in a site within the well is one. The energy 
of the two sites, the barrier and one site of the well with 
one boson, is A. This energy can be compensated if the 
barrier particle jumps to the single occupied site, so the 
energy of the two sites is U, and the transition happens 
for A ~ t7. Note that other possibilities are more expen¬ 
sive, which is corroborated by the phase diagram of the 
AB 2 chain shown in Fig. sni where the gap tends to sat¬ 
urate for larger values of A, banning quantum transitions 
in this region. A quantum phase transition around U for 
non-integer densities is a newly discovered fact, due to the 
AB 2 chain. 

In Figures 0c) and0d), we see that for the density 
p = 5/3 the arrangement of the particles passes from one 
particle in the barrier and four in the wells to five particles 
occupying two sites in the wells, and the barrier site is 
empty when X/t varies from 10 to 24, respectively. This 
situation is similar to the AB chain case for the density 
p = 3/2; therefore the quantum phase transition takes 
place around 2U. 

Finally, for X/t = 15 and p = 2, the configuration for 
the first CDW phase is shown in Fig. 0e), with five parti¬ 
cles in the wells and one in the barriers, but at X/t = 24, 
Fig. 0f) shows six particles located at two sites with the 
barrier site empty. The values of A indicate that the quan¬ 
tum transition happens around 2t7, a value for which there 
is no quantum transition in the AB chain case. On the ba¬ 
sis of Fig. 0e), note that the smallest number of particles 
in a site of the wells is two, and the barrier is always occu- 
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Fig. 10. (color online). Phase diagram for bosons in a su¬ 
perlattice type AB 2 at the thermodynamic limit. The points 
(open circles) represent the boundaries in the density matrix 
renormalization group method. 


pied. The energy for two sites, the barrier and one site of 
the well with two bosons, is U + X, while 3U is the energy 
of the two sites after the boson at the barrier jumps to the 
site with two bosons. The fluctuations will increase when 
A Ri 2U, and a superfluid region appears between the two 
CDW phases. The above discussion allows us to conclude 
that for a ABjv-i chain with a integer density p, there 
are p+1 insulator phases separated by superfluid regions; 
p insulator phases are CDW, and one is Mott insulator. 
The values of A around which the transitions occur will 
be determined by the values of N and p. 

We show the phase diagram at the thermodynamic 
limit for bosons in AB 2 chain in Fig. 1101 For commen¬ 
surate densities lower than one, the system exhibits two 
insulator phases separated by a superfluid one; the size 
of these CDW regions increases as a function of A, but 
quickly saturates. The critical values at which the CDW 
phase appears for the densities p = 1/3 and p = 2/3 are 
X/t = 0.5 and X/t = 0.1, respectively. Note that these 
phases remain stable for all values of A considered in this 
paper; i. e. for these densities the arrangement of the par¬ 
ticles is always the same regardless of the depth of the 
potential. 

The phase diagram shows us that for non-integer den¬ 
sities larger than one, there are two CDW phases, similar 
to the AB chain case, but for the AB 2 chain, we obtained 
that the transition regions depend on the density. How¬ 
ever, the transitions happen around multiples of U. When 
A = 0, the ground state is superfluid for the densities 
p = 4/3 and p = 5/3. A quantum phase transition to a 
CDW phase takes place as A increases and we found that 
the critical values are X/t = 0.1 and X/t = 0.3 for the 
densities p = 4/3 and p = 5/3, respectively. 

For integer densities, we observe that the ground state 
is Mott insulator type at A = 0 for any density p. When A 
increases, we see that CDW phases can appear, specifically 
p CDW phases for a fixed global density p (see Fig. [TUI). 




Fig. 11. Average entropy versus A for a AB chain with density 
p — 3/2 (a), AB 2 chain with p = 4/3 (b), and AB 2 chain with 
p = 5/3 (c). 

In the present paper, we consider bosons in an inho¬ 
mogeneous lattice. Information about the system can be 
obtained using the tools of information theory, which have 
been shown to be useful for determining critical points of 
systems with a harmonic potential, showing that there is 
a one-to-one correspondence between the local von Neu¬ 
mann entropy and the charge fluctuations m- Franga and 
Capelle m showed that the average of local von Neu¬ 
mann entropy is able to indicate the critical points where 
the phase changes happen in an inhomogeneous system, 
which is given by: 

i 

where = ~Trcri log 2 Ui is the local von Neumann 

entropy, and Ui = TrB<J is the density matrix of a single 
site located at i, where B represents the environment with 
L — 1 sites, and a is the density matrix of the whole system. 

In Fig. 111! we show the average von Neumann entropy 
versus X/t. In (a) we consider the AB chain, and the AB 2 
chain in (b) and (c). In all cases, we consider densities be¬ 
tween 1 and 2, and we observe the same behavior: the von 
Neumann entropy increases, reaching a maximum value, 
and then decreases. We note that it is impossible to de¬ 
termine the critical points that delimit the superfluid re- 
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gions; however we can obtain information. In all cases, for 
values of X/t on the left side, the system has a characteris¬ 
tic distribution of particles, for instance {2,1,2,1,2,1.} 

for p = 3/2 on the AB chain. The average von Neumann 
entropy has a finite value associated with the number of 
degrees of freedom. When X/t increases, fluctuations also 
increase, and we expect that the average von Neumann 
entropy will grow. This tendency continues until the sys¬ 
tem reaches its maximum number of degrees of freedom, 
which would be associated with a coherent behavior of the 
particles, i. e. with a superfluid state. Note that the posi¬ 
tions of the maximum values of the average von Neumann 
entropy agree with the middle point of the superfluid re¬ 
gions in the phase diagram for each density. Further in¬ 
crease in X/t lets the particles localize, and a different 
arrangement of the particles begins. Therefore, the num¬ 
ber of degrees of freedom decreases and the average von 
Neumann entropy is smaller. An important fact is that the 
slopes before and after the maximum are always different, 
which indicates that the insulator regions surrounding the 
superfluid phase are different. For instance, for p = 3/2, 
after the superfluid region the distribution of particles is 
{3,0, 3,0, 3,0.}. This behavior was also obtained for in¬ 

teger densities. 

During the final editing of this paper, we became aware 
of a Li et al. manuscript m. who studied a boson system 
with Xi = Acos(2(i -|- l)7r/3) using quantum Monte Carlo 
simulations. They found a phase diagram which is quali¬ 
tatively similar to the AB 2 phase diagram reported here, 
i. e. they show insulator phases for commensurate den¬ 
sities, multiples of 1/3, and observe that for p = 1 and 
p = 4/3 a superfluid region separates two insulator ones. 
They analyzed the topological properties of the insulators 
and found topological nontrivial and trivial states; specif¬ 
ically, the Mott insulators are trivial. 

The phase diagrams shown in the two papers are com¬ 
plementary, and we will use the Berry phase to determine 
the boundaries of the intermediate superfluid phase for 
densities p > 1, which results will be reported soon J!M\ . 

5 Conclusions 

Using the density matrix renormalization group method, 
we determined the chemical potential of AB^-i chains 
at the thermodynamic limit and found the phase diagram 
for the N = 2 and 3 cases. For a small energy difference 
between the A and B sites (A), we observe that insulator 
regions with peculiar charge distribution (CDW) appear 
for densities p = a/N (a an integer), except at integer 
densities, for which the Mott insulator phase still appears. 
The size of the new CDW phases for densities less than 
one grows with A, but then it stabilizes, and these phases 
remain in the phase diagram. On the other hand, for den¬ 
sities p > 1, we always found superfluid regions that sep¬ 
arate two insulator phases, a result that was confirmed 
using the von Neumann entropy. Regardless of the value 
of N, we found that for integer densities p, there are p -f 1 
insulator phases, these being p CDW phases and one Mott 
insulator phase with p particles per site. 


For ABpf-i chains, we observed that there are two 
CDW phases for any non-integer densities larger than one. 
This result can be generalized by saying that for non¬ 
integer densities larger than one, there are p -f 1 CDW 
phases, where p corresponds to the previous integer den¬ 
sity. For instance, for commensurate densities larger than 
two, we expect three CDW phases. 

The superfluid regions that separate the CDW phases 
occur around multiples of the local repulsion U, but the 
specific values and the size of these regions depend on the 
global density p and N. 
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